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ABSTRACT 

We present the results of high spatial resolution observations at 1.1 mm waveband, 
with the Submillimeter Array (SMA), towards the protocluster G10.6-0.4. The 1.1 mm 
continuum emission reveals seven dense cores, in which infall motions are all detected 
with the red-shifted absorption dips in HCN (3-2) line. This is the first time that 
infall is seen towards multiple sources in a protocluster. The dense core located in the 
center has the largest mass and mass infall rate and dominates the gas accretion of 
the protocluster, while the off-center cores have much smaller masses and mass infall 
rates, which favors the "competitive accretion" model. We also identified four infrared 
point sources in this region, which are most likely Class 0/1 protostars. Two jet-like 
structures are also identified from Spitzcr/IRAC image. 
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1 INTRODUCTION 

There are two most promising models to account for high- 
mass star formation (M >8 Mq). One is called the "mono- 
lithic collapse and di sk accretion" (|York fc Sonnhalterll2002l ; 
iMcKee fc Tar] l2003lh a nd the other is "competitive ac- 
cretion" (| Bonn efi et all 120021 : iBonnell. Vine fc Bate! 12004 
iBonnell fc Bate! 2006). The former one suggests that high- 
mass stars form directly from isolated massive gas clumps 
like low-mass star formation but with a much larger accre- 
tion rate. The latter claims that high-mass stars form in the 
center of a cluster through competitive gathering gas mass 
with the other off-center protostars. 

The competitive accretion model was proposed based 
on the observational fact that most stars are formed in clus- 
ters, and m ost massive stars re side in the center, i.e, mass 
segregation |Bonnell et al Jl2002h ■ In this model, the mutual 
potential of the protostars in the cluster funnels gas down to 
the center, which results in higher gas densities and higher 
accretion rates for the few protostars located there. Thus the 
protostars located in the center can accrete more and more 
gas and form massive stars. However, till now, no convincing 
observational evidences support the "competitive accretion" 
model. What is the mass segregation in protoclusters? How 
do the protostars interact with each other? Do the protostars 
competitively accumulate masses under a common poten- 
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tial? To distinguish these two models, high spatial resolution 
observations towards protoclusters are critically needed. 

Locating at a distance of 6 kpc, the dense clump 
G10.6-0.4 has been widely studied from centimeter to sub- 
millimeter wavele n gths at various spatial resolutions (e.g., 
iLiu. Zhang fc Hoi ( 201 lh and the references therein). The 
previous centimeter observations have pointed out that 
a burst of massive star formation in its central region 
llHo, Klein fc Haschickl Il986l : iKeto. Ho fc Haschickl 1 19871 . 
ll988l : ISollins et al.ll2005h . A flatten rotating torus is found in 
the c e ntral parsec regio n of the dense clump (|Ho fc Haschickl 
1 19861 : ILiu et all l2010al ) . Gas infalling was also revealed in 
the form of redshifted absorption seen against the contin- 
uum source i|Keto, Ho fc Haschickl H9871 . 1 19881 : ISollins et all 
120051 ). Multiple outfl ows were detected in th e 12 CO (2-1) 
and HCN (3-2) lines i|Liu. Ho fc Zhandl2010bT ). All the pre- 
vious observations indi cate a high-mass star clu ster is form- 
ing in a dense clump l|Liu. Zhang fc Holl201ll ). Thus it is 
a good target for testing different models for high mass 
star foramtion with sufficiently high spatial resolution ob- 
servations and sensitive molecu lar tracers for infall motion . 
In previous VLA observ ations (|Keto. Ho fc Haschickl Il987l . 
1 19881 : ISollins et al"1l2005t ), the NH 3 lines used as a tracer for 
infall motions only revealed the gas flow in the inner warm 
part near the central ultracompact Hn region, which cannot 
provide information of the infall motions in the other less 
dense and colder parts of the dense clump. In this paper, 
with the HCN (3-2) line , which is a good tra cer for infall 
motions (|Liu et all 1201 lal lbl: IWu fc Evansll2003r t, we can re- 
veal the pattern of infall motions in the whole dense clump 
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and uncover the accretion mode for the formation of a high- 
mass star cluster. 



2 OBSERVATIONS 

The observations of G10.6+0.4 were carried out with the 
SMA in July 2005 with seven antennas in its compact 
configuration and in September 2005 with six antennas in 
its extended configuration. The phase reference center of 
both observations was R.A.(J2000) = 18 h f0 m 28.698 s and 
DEC.(J2000) = -19°55'48.68". The 345 GHz receivers were 
tuned to 265 GHz for the lower sideband (LSB) and 275 GHz 
for the upper sideband (USB). The frequency spacing across 
the spectral band is 0.8125 MHz or ~1 km s" 1 for both con- 
figurations. 

In both configuration observations, Jupiter, Uranus and 
QSO 3c454.3 were observed for antenna-based bandpass cor- 
rection. QSOs 1741-038 and 1908-201 were employed for 
antenna-based gain correction. Uranus was observed for flux- 
density calibration. 

Miriad was employed for calibration and imaging 
|Sault. Teuben fe Wright] QUI). The 1.1 mm continuum 
data were acquired by averaging all the line-free chan- 
nels over both the 2 GHz of upper and lower spectral 
bands. MIRIAD task "selfcal" was employed to perform self- 
calibration on the continuum data. The gain solutions from 
the self-calibration were applied to the line data. 

The synthesized beam size and 1 a rms of the dust 
emission observed in the compact configuration is 3". 47 x 
2". 89 (PA=-45.7°) and 30 mjy beam" 1 , respectively. The 
continuum data combined from both configurations yield a 
synthesized beam of l".24x l".ll, P.A.=-61°.6, and 1 a rms 
of 10 mjy beam -1 in the uniform weighted map. 

IRAC data were also retrieved from the database of 
GLIMPSE. 



3 RESULTS 

3.1 1.1 mm continuum emission 

The 1.1 mm continuum image obtained from the SMA com- 
pact array is shown in red solid contours in the upper 
panel of Figure 1. Despite of the central elongated clump, 
two small lumps were found in the nort h-east and north- 
west. T he two UC Hn regions identified bv lLiu. Ho fc Zhand 
(|2010bl ) were denoted as black plus symbols. The emission 
peak of the 1.1 mm continuum emission coincides with the 
central bright UC Hn region. The deconvolved size of the 
large clump is 3".67x l".07 (PA.=-64.1°). By combining the 
data of the compact and extended configurations of SMA, 
we obtained a much higher spatial resolution (l".24 x l".ll, 
P.A.=-61°.6) image of the 1.1 mm dust continuum emission. 
As shown in the lower panel of Figure 1, The large clump 
detected by the compact array is further resolved into seven 
dense cores, which are denoted from "A" to "G". The posi- 
tions, integrated intensity and radii of these cores are listed 
in columns 2-4 of Table 1. 

Assuming that the dust emission is optically thin 
and the dust temperatur e equals the rotatio nal tempera- 
ture (87 K) of CHfCN |Beltran et al.ll201ll ). the masses 
of seven cores can be obtained with the formula M = 
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Figure 1. Upper: The flux ratio map of [4.5]/[3.6] is shown in 
color scale. The white dashed contours represent the distribu- 
tion of 8 urn emission. The contour levels are from 400 to 1200 
mjy/sr in a step of 100 mjy/sr. The 1.1 mm continuum from the 
compact configuration is shown as red solid contours. The con- 
tour levels are (-5,5,10,20,40,80,160) xO.03 Jy/beam (1 a). The 
two black plus symbols mark the positions of the UC Hn region. 
The red open circles and open boxes mark the locations of the 
blucsh ifted and redshifte d 12 CO outflow interaction signatures 
jLiu. Ho fe Zhand [ioiubl'l . The red solid circle is the location of 
the HCN (3—2) outflow. The names of the four infrared point 
sources are labeled as IRS 1-4. The two dashed long lines denote 
the directions of finger-like structures. Lower: The 1.1 mm contin- 
uum from combining the compact and extended configurations is 
shown in white solid contours and in color scale. The contour lev- 
els arc (-3,3,6,9,12,18,24,36,48,60,72,84,96,120,144,168) xO. 01 Jy 
beam -1 (1 a). The positions of the cores are marked with 
"crosses" and their names are labeled from "A" to "G". 



S U D 2 /k, u B„(T,i), where S v is the flux of the dust emission 
at 1.1 mm, D is the distance, and B v (Td) is the Planck func- 
tion. The dust opacity k v = Ki3oo(— x— — )^=0.011, where 
13=1.5 is the dust opacity index and K1300 = 0.009cm 2 g -1 is 
the dust opacity at 1 300 |xm derived from a gas /d ust model 
with thin ice mantles |Ossenkopf fc Henningj 1994). Here the 
ratio of gas to dust is taken as 100. Since the central core 
(D) is a UC Hn region, which should have a higher temper- 
ature than the clump averaged temperature derived from 
CHj 3 CN lines, we also derived a lower limit for the mass 
and volume density for it by assuming a dust temperature 
of 200 K. The inferred masses and volume density of these 
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Table 1. Parameters of each dense core. 



Name 


Offset 


s„ 


R 


n 


M 


V; sr 


CTjvt 




(",") 


(Jy) 


(10 3 AU) 


(10 7 cm" 3 ) 


(M ) 


(km s _1 ) 


(km s — 1 ) 


A 


(6.0,1-3) 


0.37(0.06) 


4.3(0.5) 


1.2(0.5) 


31(5) 


-7.50(0.06) 


1.72(0.06) 


13 


(2.1,0.7) 


0.23(0.03) 


2.7(0.3) 


3.1(1.1) 


19(2) 


-3.73(0.23) 


1.80(0.20) 


C 


(2.5,-3.5) 


0.56(0.03) 


3.5(0.2) 


3.3(0.6) 


46(2) 


-5.10(0.07) 


1.43(0.07) 


D 


(-0.3,-1.2) 


4.81(0.31^ 


6.0(0.3) 


[2.5,5.9] 


[165,400] 


-2.74(0.05) 


2.22(0.04) 


E 


(-4.3,3.0) 


0.50(0.07) 


6.4(0.8) 


0.6(0.2) 


41(6) 


-0.07(0.14) 


3.13(0.13) 


F 


(-3.7,0.4) 


0.60(0.07) 


3.7(0.4) 


3.1(1.1) 


50(6) 


-1.07(0.11) 


2.29(0.11) 


G 


(-4.9,-1.0) 


0.37(0.04) 


3.5(0.3) 


2.4(0.7) 


31(3) 


-2.35(0.15) 


1.92(0.13) 



a Since the integrated flux density at 1.3 cm is as high as 2.5 Jy dSollins et al ■l2005h . the free- free contribution 
to the 1.1 mm can't be ignored. Assuming optically thin free-free emission (S v oc v~ 01 ) iBeltran et al.l201lT ), 
1.99 Jy of free-free emission was subtracted from the total flux 6.8 Jy. 



cores are listed in columns 5 and 6 of Table 1. The masses 
of the cores range from 19 to 165/400 Mq. 

3.2 Infrared emission 

Four infrared point sources locating at the central part are 
revealed by the 8 urn emission presented as white dashed 
contours in the upper panel of Figure 1, which are denoted 
as IRS 1-4. IRS1 locates west and is not associated with the 
1.1 mm dust emission. The other three infrared point sources 
form a triangular system, whose geometric center coincides 
with the 1.1 mm emission peak. All the four infrared point 
sources have colors [3.6]- [4.5] ^0.4 and [5.8]-[8.0]^l.l, indi- 
cating they are most likely protostars (Class 0/1 objects) 
ilQiu et alj|2008l ). 

The flux ratio map of [4.5] to [3.6] urn emission is shown 
in color scale in the upper panel of Figure 1. Two finger-like 
structures are revealed and marked by the two long dashed 
lines. Despite of the finger-like structures, two red spots with 
ratios larger than 4 are found near the 1.1 mm emission 
peak. The flux ratio of [4.5] to [3.6] in the jets is comparable 
or hi gher than ~1.5 in contrast to the stars (|Takami et al.l 
|2010| ). The two finger-like structures may depict the jets gen- 
erated from the protocluster G10.6-0.4. The roots of the two 
"fingers" connect with two red spots, indicating the inter- 
action between the shocked jets and the envelope. We also 
noticed that the 12 CO outflow knots distribute along the 
finger-like structures and mainly locate at the outer layer 
or the tip of the jets, indicating the 12 CO outflow gas is 
entrained and accelerated by the collimated jets. 



3.3 Line emission 

The source averaged HCN (3-2) and CH 3 OH (5 2 ,3-4i, 3 ) lines 
are shown in Figure 2. As shown in the last panel, the HCN 
(3-2) line averaged over the whole clump has an inverse P- 
cygni profile, indicating that the whole protocluster is in 
overall collapse. As marked by the red dashed lines, red- 
shifted absorption dips were detected in the HCN (3-2) line, 
indicating gas infalling is taking place towards these dense 
cores. This is the first time that infall is seen towards mul- 
tiple sources in a protocluster. Only core "F" and core "G" 
show strong high velocity blue w ings, suggesting that the 
previous discovered HCN outflow (jLiu. Ho fc Zhandl2010bl ) 
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Figure 2. Spectra averaged over each core, green: HCN (3—2), 
red: CH3OH (52,3-41,3). The last panel presents the clump av- 
eraged spectra obtained only from compact configuration. The 
black dashed lines mark the position (-3 km s -1 ) of the systemic 
velocity. The red dashed lines mark the absorption dips. The best 
"two-layer" model fits towards the dense cores are exhibited as 
blue solid lines. 



in this region is driven by the forming protostars in core "F" 
and core "G". The CH3OH (52,3-41,3) line is single-peaked 
towards all the cores. The systemic velocities of the CH3OH 
(52,3-4i,3) line are listed in column 7 of Table 1. We found 
the CH3OH (52,3-4i,3) line is red-shifted towards the west- 
ern cores ("E", "F", and "G") and blue-shifted towards the 
eastern cores ("A", "B", and "C"). 
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Table 2. Best fitting results for the HCN (3—2) line with the enhanced two layer models 



Name 


TO 


f 


J c 


J/ 




V; sr 


(j 




M a bs 


M BH 








(K) 


(K) 


(K) 


(km s- 1 ) 


(km s- 1 ) 


(km s~ 1 ) 


(10- 2 MQ-yr- 1 ) 


(10- 2 MQ-yr- 1 ) 


B 


0.5 


0.4 


49.3 


8.1 


38.7 


-3.9 


2.5 


1.5 


0.7(0.3) 


0.5(0.3) 


C 


2.8 


0.5 


54.8 


25.2 


51.0 


-2.3 


3.6 


1.6 


1.3(0.3) 


6.3(1.6) 


D 


0.5 


0.1 


91.3 


22.8 


70.8 


-2.5 


2.7 


2.8 


[4.9,11.9] 


[16.0,243.2] 


F 


0.7 


0.3 


65.9 


8.1 


63.6 


-2.5 


2.7 


2.7 


2.3(1.0) 


1.8(0.8) 


G 


1.1 


0.4 


44.6 


11.4 


32.0 


-1.6 


2.3 


2.5 


1.4(0.5) 


0.9(0.4) 



4 DISCUSSION 

4.1 Gas inflow in the protocluster G10.6-0.4 

The redshifted absorption in t he HCN (3-2) lin e indicate 
gas infalling towards the cores l| Welch et al.ll 19871 ). Interfer- 
ometers can resolve out the extended emission and create 
kind of artifacts that may appear like absorption features, 
especially close to the line center. However, the absorption 
dips in the HCN (3-2) line are redshifted by several km s~ 
from the systemic velocity, which should not be caused by 
the missing flux. Especially, the HCN (3-2) line towards core 
D, F and G shows clearly inverse P-cygni profiles with the 
absorption dips ~5 K (~ 5cr) below the continuum, indicat- 
ing gas infalling to wards the bright continuum background 
I Welch et al.ll 19871 ). 

To investigate the properties of the absorbing gas, 
the profiles of the HCN (3-2) line were mod eled with an 
enhanced version of the "two - layer" method (|Mvers et alj 
ll99fj ; |Pi Francesco et"afll2002j ). In that model, a continuum 
source is located in between the two layers. Each layer has 
peak optical depth to, velocity dispersion a, and approach 
speed Vj„. The line brigh tness temperature at velocity V is 
jDi Francesco et aljliooa l: 



AT E 



where 



and 



(J/-J cr )[l 



• exp{-Tf)] + (1 - $)(J r 
X [1 — exp(—T r - 



-Jo) 

Tf)] 



Jcr = $J C + (1 - $) Jr 

{V -V m -V LSR ) 2 , 



Tf = Toexp[ 



2(7 2 



r exp[ 



-(V + Vin-VLSR) 2 
2<7 2 



J c , Jf, J r , Jb are the Planck temperatures of the contin- 
uum source, the "front" layer, the "rear" layer and the cos- 
mic background radiation, respectively. J = 



k exp(T /T)-l 

is related to the blackbody temperature T at frequency v, 
where h is Planck's constant, and k is Boltzmann's constant 
|Mv ers et al1ll996l ). <3? is the filling factor of the continuum 
source. We did not take outflow into the modeling. The lines 
towards core A and core E were not modeled because they 
seem to be greatly affected by the missing flux. The best fits 
are displayed in Figure 2 and the parameters of the best fits 
are summarized in Table 2. 

The kinematic mass infall rate can be calculated using 
dM/ dt=4irn^Gr>iH2RinVin, where Vi„ is the infall velocity, 
pa = 1-36 and m,H 2 are the infall velocity, the mean molec- 



ular weight, and the H2 mass, respectively. The radius and 
volume density of the cores are taken as infall radius R, n 
and volume density n in the calculation. The derived mass 
infall rates are listed in the tenth column of Table 2. 

The "competitive accretion" model predicts the ac- 
cretion rates in the virializcd cluster are primarily 
determined by the Bondi-Hoyle type accretion rates 
|Bonnell et al.ll2002l ; [Bonnell fc Bate] I2OO6I ). However, this 
argument was often doubted a nd criticized by people 
jKrumholz. McKee fc Kleinl l2005h . The Bondi-Hoyle type 
accretion rates are believed too low to significantly in- 
crease the protostars' masse s due to large velocity disper - 
sion in turbulent cores (|Krumholz, McKee fc Kleinl 120051) . 
And the radiative feedback from the massive protostars 
may also prevent Bo ndi-Hoyle accretions in the cluster 
jEdear fc Clarke! 120041 ). Does the mass accretion in G10.6- 
0.4 favor the Bondi-Hoyle type accretion? The Bondi-Hoyle 
acc retion rates of the dense cores in G10.6-0.4 are calculated 
as jKrumholz. McKee fc Kleinll2005l ): 



(GM) 2 



(1) 



where G is the Gravitational constant, p and M are the 
mean density and mass of the dense cores, respectively. The 
one dimensional velocity dispersion a — \J a 2 NT + o-'j, herm . 
The one dimensional thermal velocity dispersion arherm = 
-^L^0.51 km s~\ where 



1.36 and T fe =87 K. The 



one dimensional non-thermal velocity dispersion listed in the 
last column of Table 1 is estimated from the CH3OH (52,3- 

4i, 3 ) line through a NT = ^ch 3 oh ~ mc^on ' The calcu " 
lated Bondi-Hoyle accretion rates are listed in the last col- 
umn of Table 2. One can find except core "C" and core "D", 
the Bondi-Hoyle accretion rates of the other cores are similar 
to the observed mass infall rates, indicating the accretion in 
them is Bondi-Hoyle type. The Bondi-Hoyle accretion rates 
in core "C" and core "D" are significantly larger than their 
observed mass infall rates, indicating the accretion there is 
greatly resisted by the feedback from the massive protostars 
located in the center. However, the central massive core has 
accumulated enough mass to maintain a large accretion rate. 
With the present observations, we are unable to estimate 
how many accreted gas by the dense cores falls onto the pro- 
tostars. But the central massive protostars seem to continue 
growing mass by accreting ionized gas passing throug h the 
ioniz ation boundary of the ultracompact Hn region |Ketd 
|2002|) . 
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4.2 Competitive accretion in G10.6-0.4? 

From Figure 2, one can find the absorption dips of the HCN 
(3-2) line are at nearly the same velocity, indicating the 
gas infalling towards these cores is dominated by a common 
potential. In free fall case, the effective infall radius rj n of 
the central dense core only due to gravity can be estimated 
by use of: 



2GM 
~V T ' 



(2) 



The estimated r in is ~36000 — 87000 AU, which is much 
larger than the clump size (~ 10000 AU), indicating that 
the central dense core dominates the gravity of the clump 
and thus can competitively gather gas mass with the other 
off-center protostars from the gas reservoir of the clump. 
From Table 1 and 2, one can find that the central dense 
core (core "D") has the largest mass and mass infall rate, 
but the surrounding smaller cores have much smaller mass 
and mass infall rate, indicating that the mutual potential of 
the protocluster likely funnels gas down to the center. This 
picture favors the "competitive accretion" model. 

The difference between the "competitive accretion" 
model and the "monolithic collapse and disk accretion" lies 
primarily in how and when the mass is gathered to form the 
massive star (|Bonnelll2008l ). Does it occur in the pre-stellar 
stage due to turbulence, or does it occur subsequ ently due to 
accretion in a common potential (|Bonnelll20oi ). G10. 6-0.4 
is a much evolved protocluster. But the protostars and their 
gas reservoirs (dense cores) in G10.6-0.4 are still accumulat- 
ing masses from the gas reservoir of the clump, indicating 
that the cores in G10.6-0.4 are most likely formed due to 
accretion rather than pre-existing. 

At the present time, we can not make unambivalent con- 
clusion that the protocluster in G10.6-0.4 is formed through 
competitive accretion due to uncertainties in estimating the 
masses and mass infall rates. However, we clearly revealed 
the mass segregation in this region and the cores seem to 
still accumulate masses due to a common potential. Fur- 
ther higher sensitivity and resolution observations especially 
from the ALMA can give a more thorough understanding 
of the situations in this region especially of the interaction 
among the individual cores. 
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